Charge, density and electron temperature in a molecular ultracold plasma 
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A Rydberg gas of NO entrained in a supersonic molecular beam releases electrons as it evolves to 
form an ultracold plasma. The size of this signal, compared with that extracted by the subsequent 
application of a pulsed electric field, determines the absolute magnitude of the plasma charge. This 
information, combined with the number density of ions, supports a simple thermochemical model 
that explains the evolution of the plasma to an ultracold electron temperature. 

PACS numbers: 52.55.Dy, 32.80.Ee, 33.80. Gj, 34.80.Lx 



Under suitable conditions, the quasi-neutral system of 
ions and electrons that forms an ultracold plasma can 
develop the properties of a strongly coupled fluid. This 
intriguing degree of charged-particle correlation occurs 
when the average interparticle Coulomb energy exceeds 
the thermal energy, or T > 1, where T quantifies the ratio 
of Coulomb to thermal energy in the ion and electron 
subsystem reservoirs: T^e = (q 2 /4Treoa)/kTi }e , in which 
the Wigner-Seitz radius, a, relates to the particle density 
by, §7ra 3 = l/p. 

For r 1, a plasma can form liquid- like correlations 
PP. Beyond a correlation parameter of about 180, simu- 
lated one-component plasmas undergo a phase transition 
to form body-centered-cubic crystals 0. The creation 
of an experimental plasma with properties such as these 
would provide an important laboratory analog for the 
strongly coupled systems that appear under conditions 
of extreme temperature and density, such as in quark- 
gluon plasmas, at the centres of heavy stars and in the 
plasmas generated by the interaction of high-power lasers 
with solid targets. 

Most often, ultracold plasma experiments start with 
a magneto-optical trap (MOT) in which an atomic gas 
sample with an initial temperature lower than 100 fiK 
absorbs light from a high-resolution laser to release elec- 
trons with an average energy of 10 K or less [3]. In 
the plasma that forms, three-body electron-ion recombi- 
nation drives the initial electron temperature to a value 
above the strongly coupled regime, typically T e > 30 K. 

In a novel approach, we have excited nitric oxide 
cooled in a seeded supersonic beam to form a durable 
molecular ultracold plasma with an ion density exceed- 
ing 1 x 10 12 cm -3 [4], and an electron temperature, es- 
timated from a Vlasov analysis, to fall with expansion 
from an initial value as low as 7 K to 1 K or less over a 
lifetime that extends beyond 100 /is [5] . 

Such conditions of temperature and density suggest an 
electron coupling parameter substantially greater than 1. 
If confirmed, this degree of correlation would represent a 
significant milestone in ultracold plasma research. The 
molecular beam method offers a versatile means to study 
plasmas built both on atomic and a great many molec- 
ular precursors, and the well defined propagation of the 



plasma volume in the laboratory frame presents signifi- 
cant advantages for imaging. 

The present letter reports a measurement of the ab- 
solute excess charge of this plasma, and uses this infor- 
mation to construct a simple thermochemical model for 
its evolution from a Rydberg gas. This model serves to 
confirm our estimate of the total charged-particle den- 
sity and suggests a mechanism by which this system at- 
tains such a low electron temperature. Calculations re- 
late this feature to the high plasma densities achievable 
in a molecular beam. 

Our experimental apparatus expands nitric oxide, 
seeded at 10% in helium, through a 0.5 mm diameter 
nozzle. A 1 mm diameter skimmer selects the core of 
this jet to form a molecular beam in the second stage 
of a differentially pumped electron spectrometer. Two 
frequency-doubled pulsed dye lasers intersect this beam 
8 cm beyond the skimmer, between a pair of field grids, 
labeled Gi and G2 in Fig. [TJ The first laser excites 
NO from its ground state to the v = 0, N = level of 
the A( 2 E + ) state. From this gateway state, the second 
dye laser drives the system of molecules to the n = 51, 
I = 3 Rydberg state converging to the N + — 2 rotational 
state of the ion (51/(2)), which lies well below the lowest 
v + = 0, N + = ionization limit. At sufficient density, 
these Rydberg molecules interact, leading to a sequence 
of ionization processes that form an ultracold plasma [1] . 

This laser-crossed molecular beam illumination geom- 
etry creates a 750 /im by c.a. 2 mm prolate ellipsoid ex- 
citation volume. The cross-beam density falls off rapidly 
along the long axis of this ellipsoid, and for convenience, 
we approximate the initial plasma density distribution 
as a spherical Gaussian with a radius, r$ = 375 /im. As- 
suming saturated steps of double-resonant excitation, the 
density of NO at this point in the beam yields a molec- 
ular Rydberg density as high as 5xl0 12 cm~ 3 , or a total 
of approximately 10 9 Rydberg-state NO molecules in the 
initial excitation volume. 

We observe plasma formation heralded by a small sig- 
nal of prompt electrons. A bias applied to Gi maintains 
a very weak field (-250 mV cm -1 ) between Gi and G2, 
which extracts these promptly released electrons for de- 
tection. Subsequently, the plasma moves with the veloc- 
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ity of the molecular beam toward the detector to reach 
G2, which is grounded. On passing through this plane, 
the G2 grid extracts electrons, which accelerate in the 
100 V cm -1 held between G2 and G3 to reach the de- 
tector. The resulting waveform measures the electron 
density profile of the plasma along its axis of propaga- 
tion. 

A voltage pulse applied between Gi and G 2 , before 
the plasma reaches G2, serves as a probe to gauge the 
absolute magnitude of its excess positive charge. Using 
a Behlke high-voltage switch, we elevate Gi to a prede- 
termined negative voltage for an interval of 1 /is, 6 /is 
after the lasers have fired, which sends a pulse of ex- 
tracted electrons to the detector. The grids return to 
near ground before the plasma passes through G2. For 
best reproducibility, we sometimes find it useful to pre- 
cede this pulsed-field extraction with a 5 V cm -1 reverse- 
bias pulse to sweep away electrons weakly bound to the 
low-density wings of the plasma volume. Fig. [T] shows 
a waveform displaying the prompt signal of electrons as- 
sociated with plasma formation, the electrons extracted 
from the plasma by the pulsed field applied as it travels 
between Gi and G2, and the electrons shed as the plasma 
traverses the G 2 grid. 

The prompt signal represents the Qq electrons re- 
leased by energy pooling collisions, including Rydberg- 
Rydberg Penning ionization, and subsequent electron- 
impact ionization of Rydberg molecules. Electrons lib- 
erated by these processes escape the illuminated volume 
until the potential energy of the accumulating excess pos- 
itive charge, N* y/2/n q^/Aireoro exceeds § fee^e; where 
T' e represents the temperature of these initially formed 
energetic electrons and N* is the number of positive ions 
required to reach the trapping potential threshold for a 
plasma of radius ro [7] . 

Electron-impact ionization continues in an avalanche 
[5] that cools trapped electrons [5] until this process 
exhausts their energy or the population of Rydberg 
molecules. Electrons can also gain energy by three-body 
electron-cation recombination and superelastic collisions 
that de-excite Rydberg molecules, but dissociative chan- 
nels, open for molecules, may diminish the importance of 
these latter effects. 

We see from the initial signal that prompt loss depletes 
the plasma of a certain number of electrons, producing a 
net positive space charge. It remains now to determine 
the absolute magnitude of this charge, and reconcile its 
appearance with the electron temperature evidenced by 
the rate of plasma expansion [5]. 

By the application of pulsed field, we both tip the long- 
range potential and polarize the plasma, adding a dipole 
component to the angular dependence of its potential at 
short range. A full description of these effects would re- 
quire a self-consistent solution of the Poisson equation 
for the field-distorted electron density distribution, with 
the resulting calculated potential sensitive to the selected 
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FIG. 1: (Color online) Ultracold plasma electron signal wave- 
form. Prompt pulse, reverse-bias switching transient, pulsed 
field signal (from bottom: 37.5, 50 and 100 V cm -1 ) and 
plasma electrons extracted on transmission through G2. Sim- 
ulations model pulsed-field electron signal waveforms, neglect- 
ing the baseline elevation observed at higher extraction volt- 
ages. Inset: Molecular beam path through grids to a mi- 
crochannel plate detector. 

ion distribution [10] . We neglect this complexity, and in- 
stead, for computational as well as conceptual simplicity, 
we model the potential outside the plasma as that of a 
charged, conductive sphere in an external field: 

V(r,9) = — £ rcos6> + £;o^cos6>, (1) 

where ro is the radius of the sphere, Q is the total charge 
and Eq is the externally-applied field. This approxima- 
tion in effect: i) assumes that the high charge density of 
the plasma core screens its interior electrons, and ii) col- 
lapses all the electron density beyond a critical radius to 
form a surface charge distribution from which the pulsed 
field extracts electrons. 

The plasma charge combines with the external field to 
create a saddle point that defines e cr ; t (6'), the potential 
energy an electron on the r surface must overcome to 
escape. We calculate this flux of pulsed-field extraction 
electrons in steps by considering the trapping potential 
for each element of surface area, and then integrating over 
area and electron kinetic energy for an assumed thermal 
distribution at temperature, T e . For each step, we recal- 
culate the charge and raise the lower limit, e cr it(0), in the 
integral over electron kinetic energy for each element of 
surface area. 

Elastic collisions between the trapped electrons and 
ions transfer energy inefficiently owing to their large mass 
difference, so the ions remain at a low temperature. How- 
ever, electron-electron thermalization proceeds on the 
subnanosecond timescale of a few plasma periods [TT] , 
We therefore assume that the electron kinetic energy dis- 
tribution function is always Maxwellian. Thcrmalizing 
collisions fill-in the truncation that occurs as high kinetic- 
energy electrons from the tail of the distribution leave the 
plasma. 
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To simulate the experimental data, we select an ini- 
tial plasma charge, Qo and calculate the total number of 
electrons, Q, extracted at each applied field. If the Qo 
we choose is too small, this calculation overestimates the 
pulsed field signal relative to Qo. If we choose too large a 
Qo, the predicted Q/Qo is too small. Figure [2] compares 
experimental ratios of pulsed field to prompt signals with 
simulations scaled in this way for three values of Qo- A 
least-squares fit of the simulation to the experimental 
data yields an estimated value of Qo corresponding to 
the initial loss of 9.5 x 10 5 electrons. 
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FIG. 2: Signal versus pulsed-field normalized by the plasma 
charge. Measured data (solid circles) normalized by the 
prompt signal. Calculated lines normalized by Qo chosen for 
the simulation. Line through the data points for an initial 
charge, Ni — N e , of Qo = 9.5 x 10 5 . Upper and lower lines 
assume Qo = 4.75 x 10 5 and 1.9 x 10 6 respectively. 

A sequence of two events determines this excess charge 
and the observable electron temperature of the plasma. 
Penning collisions first convert pairs of NO Rydberg 
molecules into a free NO + ion, an electron and a de- 
excited NO molecule, which very likely predissociates to 
form translationally excited neutral atoms. The excess 
energy of ionization creates an initial population of elec- 
trons with an average energy characterized by some tem- 
perature T' e , This energy enables electrons to escape un- 
til the attractive potential of the developing NO + space 
charge traps the electrons that remain. A cascade of col- 
lisions then excite and de-excite remaining NO Rydberg 
molecules, giving rise to an avalanche of electron-impact 
ionization events. Because, on balance, these electron lib- 
erating collisions consume energy, the trapped electrons 
cool and remain confined [HI H2] • 

A few assumptions enable us to develop a simple model 
for this evolution to a plasma. Without loss of gener- 
ality, we can regard the stages of Penning and electron- 
impact ionization to occur in sequence or overlap in time. 
The initial escape of Penning electrons evaporatively low- 
ers Tg, but we assume the loss of comparatively few 
confines the rest, making this effect negligible. During 
avalanche ionization, exothermic electron-Rydberg de- 
activating collisions can occur, but this energy pooling 
process simply raises the phenomenological temperature 
we ascribe to the Penning electrons. Unlike atoms, NO 
molecules in Rydberg states cascading to lower principal 



quantum number rapidly predissociate, and we assume 
that product N and O atoms leave the interaction vol- 
ume without affecting the electron or ion temperatures. 

Experimentally, Killian et al. B found that they could 
describe the fraction of electrons confined by various plas- 
mas of defined excess energy as a uniform function of the 
number of ions scaled by the number, N* , required 
at threshold in each case to trap electrons. This function, 
confirmed by simulations [7J, relates N* to N[ — N' e , the 
excess charge or number of electrons lost, and Nt, the 
number of ions (as produced in the current instance by 
Penning ionization): 

(N- - iV') 2 
N = N ( ■ ( 2 ) 

As defined above, the threshold number N* increases for 
plasmas of higher electron temperature and larger tq ra- 
dius as, 

Our experiment measures the final excess charge, Ni — 
N e , which we equate to the early time excess charge, 
N£ — N' e> on the assumption that no avalanche electrons 
escape. We can accurately estimate the total number 
of Rydberg molecules, produced at saturation, N a , from 
the properties of our molecular beam. We have previ- 
ously measured the electron temperature of the evolved 
plasma, T e by determining its expansion rate [5]. 

Starting with Ni — N e as currently measured, together 
with known N and T e , we can examine the energy bal- 
ance for various trial values of fp, the fraction of Ry- 
dberg excited NO molecules consumed by Penning ion- 
ization. Each trial value of fp sets a balance between 
N^ the number of initially formed ions, calculated as, 
N'i = fpN /2, and N", the number of excited NO 
molecules remaining to be ionized by the electron-impact 
avalanche, as 

N? = N' e ' = (l-f P )N . (4) 

Having N[, we obtain N* from Eq. [2] We can then use 
Eq. [3] to determine a temperature T' e of the Penning 
electrons. 

Each ion produced by electron impact ionization of a 
nitric oxide molecule in the 51/(2) Rydberg state con- 
sumes the binding energy AE n = 8.38 x 1CT 22 J. Thus, 
avalanche ionization to form N" ions lowers the electron 
temperature of the plasma to a value, T e , such that: 

_ \k B T e N' e - N"AE n 

lk B {N> e + N») ■ [b) 

N' e approximately equates to N[ over a wide range of 
fp, and we can use Eqs. |2]-|4]to write a simple expression 
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for T e in terms of the measured charge, N — N e , and 
the fraction, fp, of Rydberg NO molecules consumed by 
Penning ionization. 



T, 



g^Nj - N e ) 2 /V2nr 2ne - (1 - f P )N AE n 
§fc B JV (l-/p/2) 



(6) 

Fig. [3] shows values of T e predicted by Eq. [6] as a 
function of fp for assumed values of Nq, decreasing from 
a saturated N = 1.1 x 10 9 (density, ji e = 5x 10 12 cm -1 ). 




FIG. 3: Electron temperature upon evolution to plasma pre- 
dicted from thermochemical energy balance as a function of 
fp, for assumed values of initial Rydberg density, No, reading 
from the left, of 1, 2, 3, 4 and 5 x 10 12 cm -3 . Below f p = 0.05, 
Eq. [6] necessarily incorporates additional terms to account 
for the difference between N' and NL. 



This variation in plasma electron temperature reflects 
the simple thermochemical balance in which electron 
heating, manifested in the observed excess charge, dissi- 
pates in the cooling necessary to complete avalanche ion- 
ization. For the purposes of illustration, consider an exci- 
tation volume containing a population of NO molecules 
excited to the 51/(2) Rydberg state with a density of 
2 x 10 12 cm~ 3 . Under these circumstances the Wigner- 
Seitz radius is 500 nm, compared with an orbital radius 
of 130 nm. The corresponding Erlang distribution func- 
tion of nearest-neighbour distances predicts that slightly 
more than 6 percent of this population will exist as pairs 
spaced by 150 nm or less [12] • R seems reasonable to 
assume that this population of proximal molecules will 
interact to produce Penning electrons on a very short 
timescale, and that these electrons will both in part es- 
cape to create a trapping potential and initiate an en- 
dothermic electron avalanche. 

Thus, for N = 4.4 x 10 8 , an assumed fp = 0.06 yields 
N[ = 1.3 x 10 7 , which produces the measured excess 
charge, (N t -N e ) = 9.5 x 10 5 for a 7^=1600 K. Electron- 
impact ionization of the Rydberg molecules that remain 
consume 3.5 x 10~ 13 J of electron kinetic energy, lower- 
ing T e to 7 K. A point plotted on Fig. [3] indicates this 
No, fp intersection. The selected value of Nq reflects ei- 
ther a smaller than saturated initial preparation or just 
the number of NO molecules ionized when Rydberg pre- 
dissociation overtakes ionization and cooling stops. 



This simple model neglects some obvious aspects of 
the plasma evolution dynamics. For example, as men- 
tioned above it considers neither the reduction in T' e that 
occurs as a consequence of early electron escape, nor in- 
crease in T e by electron-impact deactivation of NO Ry- 
dberg states. However, we know that the plasma loses 
comparatively few hot electrons to form its trapping po- 
tential, and we can assume that collisional deactivation 
of Rydberg NO accelerates predissociation, which rapidly 
depletes the plasma of molecules in quantum states below 
the bottleneck for collisional ionization P3J [T3] . 

Regardless of its simplifying assumptions, this crude 
energy balance does serve to show that a feasible route 
exists for the molecular plasma of NO + ions and elec- 
trons formed in a molecular beam to reach a low electron 
temperature. The variation in the cooling range with Nq 
shows clearly how high density plays a critical role in 
providing a sufficient reservoir of Rydberg bound states 
to cool the system as it evolves. Fortunately, typical con- 
ditions attained in pulsed dye laser crossed differentially 
pumped molecular beams yield excited state densities in 
this range for many molecular systems. 
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